The Physics of Supernova Remnant Blast Waves. II. Electron-Ion 
Equilibration in DEM L71 in the Large Magellanic Cloud 
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ABSTRACT 

We present analysis and modeling of X-ray spectra from the blast wave shock of DEM L71 in 
the Large Magellanic Cloud. This remnant exhibits widespread Balmer-dominated emission 
characteristic of nonradiative shocks in partially neutral gas. We have used both Chandra ACIS- 
S data and optical Fabry-Perot spectra of the blast wave to measure the electron and proton 
temperatures, respectively. In principle, when combined, these measurements can determine the 
degree of electron-ion temperature equilibration (go = T e /T p ) immediately behind the shock 
front. In our X-ray analysis we fit Chandra spectra of three nested regions behind the blast wave 
under three different scenarios: (1) a planar, initially unequilibrated shock (go — m c /m p ), where 
the downstream electron and proton temperatures equilibrate through Coulomb collisions, (2) a 
planar, immediately equilibrated shock (go = 1) and (3) a spherical, equilibrated shock under 
Sedov evolution. Using independent measurements of T e and T p we find that the X-ray spectra 
from the fastest blast wave locations (V s ~ 700 — 1000 km s _1 ) are consistent with little or no 
equilibration at the shock front and are inconsistent with full equilibration. In contrast, spectra 
from regions showing slower blast wave speeds (V s ~ 400—600 km s _1 ) allow full equilibration but 
exclude zero equilibration. In order to directly constrain the initial equilibration, we incorporated 
our knowledge of the proton temperatures into our X-ray models to build planar shock models that 
allowed for a variable go. This model confirmed and strengthened the above results. Specifically, 
we found that X-ray spectra from an intermediate velocity shock (V s ~ 800 km s _1 ) were 
consistent with intermediate equilibration, excluding both go — m e /m p and go = 1 at greater 
than la. Overall, our results support the picture of decreasing electron-ion equilibration with 
increasing shock speed found from previous studies of optical spectra in other Balmer-dominated 
supernova remnants. 
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1. Introduction 

The outer blast wave of DEM L71 provides 
a valuable laboratory to address the question of 
electron-ion equilibration at collisionless shocks. 
These shocks have been traditionally explained 
by their Rankine-Hugoniot jump conditions (e.g., 
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Spitzer 1978), which quantify the partition of en- 
ergy into thermal and bulk kinetic forms. How- 
ever, for collisionless shocks the actual heating is 
not straightforward because at extremely low den- 
sities, the shock front is thinner than the mean free 
path for collisions. Therefore, the heating of parti- 
cles at the shock-front must be produced by some- 
thing else, generally assumed to be plasma waves 
(for a review see for instance, Draine & McKee 
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1993; Laming 2000). Since the heating is not col- 
lisional there is no guarantee that it will produce a 
thermal distribution, and in fact both simulations 
and in-situ observations of solar-wind shocks sug- 
gest that it does not. The plasma heating mech- 
anisms, to both thermal and cosmic-ray energies, 
are as yet poorly understood. These collisionlcss 
shocks are ubiquitous in astrophysical situations, 
but are only spatially resolved in interplanetary 
shocks (at lower speeds) and SNRs. The parti- 
tion of energy into thermal and non-thermal pop- 
ulations at collisionlcss shocks has broad reaching 
implications for the dynamics of the ISM, shocks 
from merging galaxies, and hot gas in clusters of 
galaxies. Our previous work on E0102.2— 7219 
studied the partition of energy into thermal elec- 
trons versus cosmic rays (Hughes, Rakowski, & 
Decourchelle 2000). Here we examine the parti- 
tion of energy between the two main thermal pop- 
ulations, electrons and protons. 

In a small sub-sample of supernova remnants 
(SNRs), the X-ray and radio emission from 
the blast wave is accompanied by optical (Ha, 
Balmer) line emission. The line profiles of these 
Balmer-dominated spectra provide a direct mea- 
surement of the proton temperature at the shock 
front, thus also placing limits on the shock ve- 
locity. The optical spectra of Balmer-dominated 
SNRs are produced by collisional excitation when 
partially neutral interstellar hydrogen is overrun 
by the blast wave (Chevalier & Raymond 1978). 
Like other non-radiative shocks, post-shock cool- 
ing losses are negligible for Balmer-dominated 
remnants. Each Balmer emission line consists of a 
narrow and a broad velocity component (Cheva- 
lier, Kirshncr & Raymond 1980). The narrow 
component is produced when cold H I, overrun by 
the shock, is collisionally excited by electrons and 
protons before being ionized. The cold neutral 
hydrogen can also charge exchange with fast post- 
shock protons, producing fast neutrals. These 
can then be collisionally excited, producing the 
broad velocity component (Chevalier et al. 1980). 
The Balmer emission arises in a thin (< 10 15 cm) 
ionization zone; therefore, the width of a broad 
Balmer line yields the proton temperature imme- 
diately behind the shock and hence also limits 
the range of possible shock velocities (Chevalier 
et al. 1980, Smith ct al. 1991). Combining this 
information with the broad-to-narrow flux ratio 



allows us to estimate the degree of electron-ion 
temperature equilibration, go = (T c /T p ) , due to 
plasma processes at the shock front (Tycho, SN 
1006 & RCW 86; Ghavamian et al. 2001). 

A more straightforward method of estimating 
the initial temperature equilibration in a colli- 
sionlcss shock is to combine the proton temper- 
ature estimated from the FWHM of the broad 
Balmer line with the electron temperature mea- 
sured from X-ray spectra. For this purpose, the 
Large Magellanic Cloud remnant DEM L71 is 
particularly appropriate because it is completely 
encircled by Balmer-dominated filaments (Tuohy 
ct al. 1982). We have utilized Fabry-Perot imag- 
ing spectroscopy to measure the width of the 
Ha broad component along most of the rim of 
DEM L71 (Ghavamian et al. 2003; hereafter 
GRHW03). Here we report results comparing 
T e from our Chandra analysis with T p from our 
Fabry-Perot analysis to infer the value of .go 
for multiple apertures around the blast wave of 
DEM L71 . 

Davics, Elliot, & Meaburn (1976) first identi- 
fied DEM L71 as a SNR candidate in an optical 
survey of the Large Magellanic Cloud (LMC) be- 
cause of its shell-like morphology that extended 
83" by 60" in diameter. The detection of X-ray 
emission from DEM L71 in the Einstein survey of 
the LMC by Long, Hclfand and Grabelsky (1981), 
confirmed its nature as a SNR. DEM L71 is ex- 
tremely faint in the radio band, where it has only 
been detected at 843 GHz with a flux of 10 mJy 
by Mills et al. (1984). In optical follow-up, Tuohy 
et al. (1982) detected filamentary Ha emission but 
only a few faint knots of [O III] and [S II] , and thus 
categorized it Balmer-dominated SNR. Smith 
ct al. (1991) included DEM L71 in their study 
of 6 Balmer-dominated remnants. Using longslit 
spectroscopy of the brightest portions of the shock 
in DEM L71 they estimated shock velocities of 
300—800 km s _1 from the Ha broad component 
widths. In a recent ASCA study of LMC SNRs, 
Hughes, Hayashi & Koyama (1998) fit the X-ray 
spectrum of DEM L71 and showed that the data 
were well described by a non-equilibrium ioniza- 
tion (NEI) Sedov model, particularly for the case 
of only minimal initial equilibration between elec- 
trons and ions. 

Further investigations into the nature of DEM L71 
utilizing the Chandra observation are given in two 
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companion papers. GRHW03 discusses the de- 
tails of the Fabry- Perot analysis of the blast wave, 
and implications for the evolutionary state of the 
remnant. In Hughes et al. (2003) we investigated 
general properties of the remnant and the compo- 
sition of the ejecta. DEM L71 exhibits a promi- 
nent double shock morphology with a clear spatial 
separation between the reverse shock (ejecta) and 
blast wave (ISM). The mass (~ 1.5M Q ) and com- 
position of the ejecta (Fe and Si- rich), strengthen 
the case for a Type la SN origin for DEM L71 . 

2. Observations and Data Reduction 

2.1. RFP optical observations 

Optical spectra of DEM L71 were obtained us- 
ing the Rutgers Fabry-Perot (RFP) imaging spec- 
trometer as described in GRHW03. For complete- 
ness we summarize the relevant analysis here. The 
14 RFP scans, centered on the wavelength of Ha 
at the systemic velocity of the LMC, provide essen- 
tially monochromatic images of the SNR. Ha line 
emission spectra can then be extracted from any 
position on the remnant. We defined 16 apertures 
around the blast wave, which contained approxi- 
mately equal numbers of counts and avoided any 
stars or radiative filaments. Fits of the Ha line 
profile determined the broad component widths 
and broad-to-narrow flux ratios for each aperture. 
For the purpose of comparison to the X-ray re- 
sults, larger regions were required to obtain rea- 
sonable statistics in the Chandra X-ray spectra. 
Adjacent regions with consistent broad component 
widths and broad-to-narrow ratios were combined 
to provide five final optical apertures (Figure 1, 
Table 1). The Ha line profiles from these aper- 
tures were then refitted. The proton temperature, 
T p , and shock velocity, V s , were calculated from 
the measured broad component width for the full 
range of g (GRHW03). The values computed 
for the extrema of go = m c /m p and g —I are 
reprinted here in Table 2. 

2.2. Chandra X-ray observations 

We observed DEM L71 using the back-side- 
illuminated chip (S3) of the Advanced CCD Imag- 
ing Spectrometer (ACIS-S) in full-frame timed ex- 
posure mode starting on 2000 January 04 for 45.4 
ks (OBSID 775). At this time, the ACIS fo- 
cal plane temperature was —110° C. The starting 



point of the reduction was the events lists from the 
revision 2 level of the standard processing. The 
CCD events were corrected for charge transfer in- 
efficiency (CTI) using software developed at Pcnn 
State University (Townsley et al. 2000). Once 
these corrections are applied the entire chip can 
be modeled with a single response function. Sub- 
sequently, the event file was screened for grade 
(retaining only values 0, 2-4, or 6), bad pixels, 
high background times and times with incorrect 
aspect. The final filtered events file corresponds 
to a livetime-corrected exposure of 34.5 ks. 

Using the optical image of DEM L71 as a tem- 
plate, we verified that the Chandra X-ray positions 
are accurate to ~0.5". This method was used be- 
cause none of the 16 X-ray point sources detected 
on the S3 chip were coincident with an optical 
counterpart from the USNO-A2.0 star catalog. 

The spectrum of the blast wave is dominated 
by emission below 0.8 keV. Therefore, this soft 
band was used to define the blast wave regions 
matched to each optical aperture. We adaptively 
smoothed the soft band image using IMSMOOTH 
from the Chandra Interactive Analysis of Obser- 
vations (CIAO) package (with a minimum signal- 
to-noise ratio of 4), then chose a contour level of 
0.8 counts per pixel to define the edge of the blast 
wave. Three polygonal regions, each four pixels 
wide (~ 2"), were defined inward of the outer con- 
tour towards the local center of curvature (Figure 
2). Blast wave spectra were extracted from the 
filtered events file using these nested regions. The 
three succeeding regions were constructed in order 
to determine not only the blast wave electron tem- 
perature but also its downstream evolution. The 
largest temperature gradient is expected close to 
the current shock front. This and the fact that 
shortly behind our defined regions other shock 
fronts and the interior (reverse shock) emission 
become evident, limited the distance over which 
we could extract blast wave spectra. In turn, 
this limited the number of nested regions to three 
because we required the extraction regions to be 
wider than the Chandra point-spread-function as 
well as wide enough to obtain sufficient statistics. 
Note that the length of each X-ray extraction re- 
gion was defined to match the full extent of the 
corresponding optical apertures. 

Separate ancillary response files (ARFs) for 
each spectral extraction region around the blast 
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wave were generated using standard CIAO proce- 
dures. Although these were retained in the final 
fits, we note here that no significant differences be- 
tween the ARFs were found. For all spectra, back- 
ground was taken from a large annulus outside the 
remnant covering radii l'-2'. We did not attempt 
to account for the recently discovered degradation 
in the low energy quantum efficiency of the CCD. 2 
Instead we compensate for this effect by allowing 
the column density in our fits to be a free param- 
eter. 

3. X-ray modeling 
3.1. Theory 

To measure the electron temperature at the 
blast wave, one has to make some assumptions 
about the evolution of both the temperature and 
ionization state downstream. These assumptions 
may bias the measured electron temperature and 
hence the final derived value go- To address this 
concern we chose to investigate three very differ- 
ent nonequilibrium ionization (NEI) evolutionary 
models for the blast wave: (1) a fully equilibrated 
(go — 1), constant velocity, planar shock model, 
(2) an initially unequilibratcd (go = m c /m p ) 
planar shock model with post shock equilibration 
on Coulomb collisional timescales (Spitzer 1978), 
and (3) a fully equilibrated (g — 1) Sedov model. 
In all three models there are only two important 
dynamical parameters, the current average tem- 
perature at the shock front, T ave ^ s , and the final 
ionization timescale, n c t of the innermost shock 
region. All other temperatures and timescales 
are constrained to follow the appropriate evolu- 
tion. The temperatures within a parcel of gas that 
may evolve function of time are the electron 
temperature T c , the proton temperature T p , and 
their number-density-weighted average tempera- 
ture, T ave . Note that throughout this paper, "av- 
erage temperature" is used to denote the number- 
density- weighted average over the particles, and 
not a spatial or temporal average. All ions were 
considered to have the same temperature as the 
protons for the sake of the average temperature, 
a small effect for the metal-poor LMC interstellar 
medium. 



2 see the calibration area of the Chandra 
website for the announcement and details: 
http : //exc . harvard . edu/cal/Acis/Cal_prods/qeDeg/index . 



Planar shocks are a reasonable approximation 
for (1) straight or small shock segments where 
the curvature is unimportant for the density pro- 
file, and (2) time periods that are short enough 
that the deceleration of the blast wave due to 
the accumulation of swept up material is negli- 
gible. The planar model has a constant velocity 
and hence a constant average temperature behind 
the shock. Initially, we chose to model the two 
extreme cases for the initial degree of electron- 
ion equilibration, go = 1 and go = m c /m p . For 
.go = 1, all temperatures (T e , T ave ^ s , and T p ) are 
equal and constant throughout the shock and only 
the ionization timescale varies behind the shock. 
For g = m c /m p , T c and T p slowly equilibrate to 
T aV e,s via Coulomb collisions (Spitzer 1978). We 
use the expressions for the time variation of tem- 
perature from Itoh (1978) and Cox & Anderson 
(1982) to describe the variation of T e as a function 
of n c t behind the shock. Note that while this for- 
mulation was derived for Sedov evolution, it is ac- 
tually equally appropriate for planar shocks. For 
examples of the evolution of planar shocks see sec- 
tion 5. 

Considering the advanced age of DEM L71 
(^4400 years, GRHW03) it was important to in- 
clude a Sedov solution model. Even under full 
temperature equilibration, the Sedov solution im- 
plies a spatially varying electron temperature be- 
hind the shock, which could mimic the temper- 
ature variation from Coulomb equilibration in a 
.go = m c /m p planar model. The Sedov solution 
used here assumes a spherical geometry and uni- 
form preshock medium. For simplicity, we only 
modeled the case of full electron-ion equilibra- 
tion, i.e., T c — T ave — T p . There are two mech- 
anisms that determine the evolution of T ave and 
n c t. First, the shock decelerates as it sweeps up in- 
terstellar gas, causing T ave ^ s to decrease with time. 
Second, the volume adiabatically expands as the 
shocked material moves out at (3/4)V s , i.e., the in- 
terior density and temperature decrease. Parcels 
of gas that are shocked at different radii, rj, and 
times ti started off with different T ave ^ s and V s . 
Hence each parcel will follow a unique tempera- 
ture and density evolution. 

The beauty of the Sedov solution lies in its 
self-similarity. All shock quantities can be related 
by dimensional analysis, and the interior profiles 
can be expressed in terms of dimensionless re- 



4 



duced variables. Hence it is possible to construct 
a generic set of reduced temperature and density 
histories as a function of current fractional radius 
(see Figures 3 and 4). We start with the Sedov 
solutions for the shock quantities and the interior 
profiles, herein denoted as the shock and interior 
solutions respectively (see Sedov 1959). The shock 
solutions can be written to relate the current ra- 
dius (R s ), temperature (T avejS ), velocity (V s ), and 
the swept up mass (Msu) to the explosion en- 
ergy, the hydrogen number density of the ambient 
medium (no), and the time since the explosion (t). 
Following Cox & Franco (1981), we express the in- 
terior solutions in terms of the dimensionless quan- 
tity (3 which ranges from 0.8 to 1 from the center to 
the edge of the blast wave and finely samples the 
outermost radii. 3 This variable is also convenient 
because it only depends on nondimensional com- 
binations of the parameters in the problem. These 
interior solutions describe the profiles of reduced 
variables, the ratios of physical quantities relative 
to the immediate post-shock values: radius (r), 
density (x), pressure (y), initial radius (r^), and 
the mass enclosed (fi) within that (3. The reduced 
temperature is just the ratio of the reduced pres- 
sure to the reduced density. For clarity, we note 
that for a given parcel of gas there are three times 
of importance in the formulation of the tempera- 
ture and density histories: the time at which the 
parcel was shocked, tf, the time at which the evo- 
lution is being evaluated, t; and the current epoch, 

t c . 

To derive the density and temperature histo- 
ries of a parcel of g function of t, we need 
to connect the shock and interior solutions. We 
make use of the fact that the mass enclosed be- 
hind the radius of a given parcel of gas is constant, 
i.e., that portions of the post-shock material never 
overtake the gas ahead of them. This allows us 
to match the equations for the swept up mass at 
the time the parcel was shocked, Msu(U) /Msu{t) 
with those for the fraction of the total swept-up 
mass that is enclosed within a given interior ra- 
dius, n(fi) = M enclosed /M su to get (i{U/t). This 
history defines where any given parcel was within 
the remnant as a function of time, as well as its 

3 There is a minor typographical error in the Cox & 
Franco (1981) definition of (i: the coefficient three tenths 
should be inverted to ten thirds, in order to match the so- 
lutions given by Sedov (1959). 



density and temperature relative to the shock val- 
ues. We can then combine this with the shock 
solutions which relate the shock quantities at t to 
the current shock quantities. Thus we obtain fi- 
nal solutions for the ratio of the evolving T ave and 
n c t of each parcel to the current shock values in 
terms of the ratio of t to t c . These generic tem- 
perature and ionization timescale histories can be 
multiplied by any chosen time, density and tem- 
perature to give specific solutions. 

There are a number of important features of the 
Sedov solution that can be seen in the generic tem- 
perature and ionization timescale histories shown 
in figures 3 and 4. The three curves shown corre- 
spond to the current fractional radii of the three 
X-ray spectral extraction regions behind the shock 
in DEM L71 . They represent the evolution of 
parcels of gas whose final radii are at the inner- 
most edges of these three regions. Note the steep 
decrease in temperature both in time and cur- 
rent radius. Material that now lies at 0.85r s was 
shocked at about one-fifth the lifetime of the rem- 
nant when the shock was almost six times as hot. 
Even the current temperatures across the three 
spatial zones vary by more than a factor of two. 
This leads to a spatial gradient in the electron 
temperature for the Sedov model that is far larger 
than the variation in temperature across regions 
for the minimally equilibrated (g = m c /m p ) pla- 
nar shock model, which is at most 40% (see Table 
5). 

The ionization timescale of the Sedov model 
also has a very different evolution than in the pla- 
nar case. Adiabatic expansion of the post shock 
gas causes both a decrease in temperature as well 
as a decrease in the density, reducing the rate of in- 
crease of n e t. The global maximum value of n e t is 
reached for parcels that are now at ~ 0.9i? s . How- 
ever, the temperature of the parcel which is cur- 
rently at 0.85r s was more than three times hotter 
than the current T ave ^ s for almost half of its evo- 
lution. Since the ionization rates depend on the 
temperature, the final n c t by itself is not a very 
accurate measurement of the ionization state pro- 
file behind the shock. This complicates a direct 
comparison of the n c t values between the planar 
and Sedov models. 
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3.2. NEI implementation 

Other authors have implemented non-equilibrium 
ionization in Sedov and planar geometries (Shull 
1982; Hamilton et al. 1983; Kaastra & Jansen 
1993; Hughes et al. 1998; Borkowski, Lyerly & 
Reynolds 2001) but for the sake of flexibility 
and the specific needs of the selected regions in 
DEM L71 , we chose to modify our own NEI code 
(Hughes & Singh 1994). In all cases we construct 
one evolving shock model that is then used to fit 
the three nested regions simultaneously. 

For both the go — m c /m p and g = 1 planar 
models, we divided the post-shock region into 200 
logarithmically spaced regions in n c t . Since the 
average shock temperature is constant in the pla- 
nar model, all parcels of gas will undergo the same 
electron temperature and ionization state history. 
At a constant shock velocity, the n c t steps repre- 
sent both the time since the material was shocked 
and the distance behind the shock. This allows us 
to calculate the emission at each timestep and sum 
them to find the total emission from the shock. For 
each step we calculated the ionization state and 
emission using the electron temperature at that 
time determined by the solution of Itoh (1978) and 
advanced the ionization state of the material dur- 
ing that n e t interval as in Hughes & Singh (1994). 
The final ionization state from each step was then 
used as the initial ionization state for the follow- 
ing interval. The emissions were summed, weight- 
ing by An c t, which for a planar shock is also the 
width perpendicular to the shock. For the purpose 
of fitting the nested DEM L71 blast wave spectra, 
each model was cut into 3 sections linearly spaced 
in n c t . We assumed that the current blast wave 
lies at the outer edge of the outermost aperture. 
However, the fitted temperatures were not partic- 
ularly sensitive to the exact choice of blast wave 
location. Finally, we included a simple spherical 
projection of the three modeled spatial regions to 
account for the contribution of the outer emission 
zones to the inner spectral regions. The projection 
in this case was done only for the three large data 
extraction regions and not the individual n c t in- 
tervals, so that the different values of R s for each 
aperture could be accounted for during the fits. 

For the Sedov NEI model, first the generic his- 
tories were scaled to a given current shock temper- 
ature and product of the ambient density with the 



remnant age. The ionization state of the material 
at every radius (parcel) was calculated separately 
using its T e and n e t history, then the emission was 
computed for the current temperature of that par- 
cel. A grid of 100 points in (3 (over its full allowed 
range of 0.8 to 1.0) was used for both the time 
evolution histories and the emission as a function 
of radius. A single set of nominal fractional radii 
(0.85, 0.89, 0.95, calculated for an 8.8 pc current 
radius) was used for all regions in DEM L71, de- 
spite their slightly different radii (see GRHW03 
figures 7 and 8), so that a common grid of models 
could be applied to all spectra. At these nominal 
fractional radii spherical projection of every indi- 
vidual segment in j3 was used to weigh the contri- 
bution of any given parcel to each of the nested 
regions behind the shock. Note that the Sedov 
model is inherently spherical, in contrast to the 
case described in the preceding paragraph where 
we were simply adding a spherical projection to 
an otherwise planar model. 

4. X-ray spectral fits 

In the preceding section, we introduced the dy- 
namical variables, T avf , tS and n c t that we wish to 
measure with our X-ray fits. However, the spec- 
tra also depend on the composition of the gas and 
the absorbing column density to the SNR. Ideally, 
one would like to allow the abundance and col- 
umn density parameters to be free for each region 
to allow for variations in the composition of the in- 
terstellar medium around the rim of the remnant, 
and any local absorption. Furthermore, assuming 
a particular abundance set may bias the temper- 
ature measurements, because they are not inde- 
pendent variables. For example, choosing lower 
abundances will often bias the fits towards higher 
temperatures. Unfortunately our preliminary fits 
revealed that the current statistics are insufficient 
to constrain these additional parameters. 

For both planar and Sedov fits the abundances 
of O, Ne and Fe were at first allowed to be free, 
with all other species fixed at the LMC abun- 
dances of Hughes et al. (1998). The absorbing 
column density, Ah, was also free in these pre- 
liminary fits. The abundances included in our 
column density parameter were kept at solar val- 
ues. Since the total column densities we found 
were similar to the minimum absorbing gas col- 
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limn in the Galaxy along the line of sight to the 
LMC (5 - 6 x 10 20 atoms cm" 2 , Hciles & Cleary 
1979), no additional LMC component was deemed 
necessary. From these preliminary fits we deter- 
mined the average abundances and column density 
across all regions and all models. Most regions did 
not have high enough statistics to constrain the 
individual abundances. Furthermore, the F-test 
showed that freeing the abundance and pa- 
rameters did not significantly improve the fit over 
using the average values. Hence a common set of 
abundances for all models and regions was used in 
order to allow a more straightforward comparison 
both between different regions and for a single re- 
gion under different models. In the final fits we 
fixed the abundances and column density to three 
sets: the average values and the average values 
plus or minus their root mean square (RMS) devi- 
ations. The final values are listed in Table 3. They 
are all reasonable for the interstellar medium in 
the LMC compared to the abundances of Hughes 
ct al. (1998). We fitted the spectra using the RMS 
deviations in the column density and abundances 
in order to incorporate the uncertainty in the tem- 
perature and timescale due to (1) the uncertainty 
in our estimate of the average abundances and col- 
umn and (2) possible variations in abundances and 
column density around the outer rim of the rem- 
nant. This estimate of the errors is quite conser- 
vative since the RMS errors are much larger than 
the errors in any given preliminary fit and include 
any variations due to the choice of model. 

The uncertainty in the abundance and column 
density parameters prevented us from distinguish- 
ing the three shock models on the basis of their 
X 2 values. Not only was every model the "best" 
choice for some region based on x 2 > but the varia- 
tion in x 2 across models was always smaller than 
across the RMS deviations in abundance and col- 
umn density (not shown). 

The nested spectra for each of the five aper- 
tures are shown in Figure 5, plotted against the 
go = m c /m p , planar model for comparison. The 
best-fit parameters of the final fits for all regions 
and all shock models are given in Table 4. The av- 
erage shock temperatures predicted by the three 
shock models are very different because the X-ray 
emission depends primarily on the current electron 
temperature. Roughly speaking then, the values 
for T aveyS are just those that give the correct T e for 



the "dominant" region in any given aperture under 
the assumptions of that shock model. The implied 
average shock temperatures are useful, however, 
for comparison with the optically determined pro- 
ton temperatures. 

Under the assumption of full equilibration, g = 
1, the electron, proton and average temperatures 
are all equal at the shock. Hence if either the 
planar or the Sedov g = 1 models are correct 
the optically measured T p should match the X- 
ray fitted T ave ^ s . This is the case only for region 
X5, the blast wave region with the slowest shock 
speed. For all other regions T p is much greater 
than the T ave ^ s values for full equilibration, im- 
plying that the electron temperature at the shock 
is lower than the proton temperature. The go = 
m e /m p planar fits confirm this result. Under the 
g = m c /m p assumption, T p ~ (n c + n p )T ave /n p 
or T p ~ (2.3/1. l)T ave . For all regions except X5, 
the X-ray implied T p is consistent with the Ha 
proton temperature within the la errors. This 
agreement occurs on the low side of our allowed 
T ave ^ s range hinting that partial equilibration may 
be preferred. 

All three NEI models indicate the same gen- 
eral result: by comparing the average shock tem- 
peratures to the Ha proton temperatures, we can 
already exclude the case of full electron-ion equi- 
libration in four out of five apertures, indepen- 
dent of the global model. In contrast, region X5, 
which has the slowest shock speed, is consistent 
with go — 1 but not with go — n^e/f^ip under the 
assumption of a constant velocity shock, and is 
marginally consistent with go = 1 under the Se- 
dov model. 

In addition to a comparison with the optically 
measured proton temperature, the parameters of 
the X-ray fit also imply a certain evolution of the 
electron temperature and ionization state. Table 
5 lists the implied values of T e and n e t for each of 
the nested regions behind the shock in the three 
NEI models. The different model assumptions 
lead to contrasting spatial gradients in the elec- 
tron temperatures. The ionization timescale val- 
ues also differ strongly from model to model, but 
this can be understood in terms of their underlying 
assumptions. For the go = m e /m p constant veloc- 
ity model the fitted n e t values are larger by factors 
of 10% to 100% than those for the go — ^ constant 
velocity model. Early on the electron tempera- 
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tures were lower and hence ionization progressed 
less rapidly. The fits therefore require a larger n c t 
value to reach the same ionization state. For the 
fully equilibrated Sedov model, the final ionization 
timescales are much lower than the go = 1 planar 
case, because earlier the material was much hotter 
than now, so the ionization occurred more rapidly. 

Deeper observations, with sufficient statistics 
to constrain the electron temperature and ioniza- 
tion state in each of the nested regions separately, 
would clearly be able to discriminate between the 
three models given here, based on their strong dif- 
ferences in the implied evolutions. 

5. Estimates of the Initial Electron-Ion 
Temperature Equilibration 

In the previous section we showed that our inde- 
pendent measurements of the electron and proton 
temperatures led to a consistent picture of the ini- 
tial temperature equilibration at each blast wave 
region regardless of which NEI model was used. 
Here we will assume a constant velocity shock and 
present an estimate of g from the independent T e 
and T p measurements, and a direct determination 
of go using the proton temperature to constrain 
our X-ray models. 

In figure 6 we plot the electron temperature as a 
function of ionization timescale for three types of 
constant velocity shock models, using region XI 
for illustration. The diamonds denote the elec- 
tron temperatures and ionization timescales im- 
plied by the fitted values of T ave ^ s and n c t for the 
fully equilibrated X-ray model, while the crosses 
are those for the go — m c /m p model. The curves 
plot the Coulomb evolution of the electron temper- 
ature given the best-fit value of the proton temper- 
ature, and the full range of possible go values. All 
curves eventually equilibrate (flatten out) to the 
average temperature. The average temperature is 
not the same for all go values, because the initial 
proton temperature is kept constant: T avf , tS will 
be lower if the electrons are not fully heated ini- 
tially. Hence, the top curve represents full initial 
temperature equilibration, where T e , T p , and T ave 
are equal and remain constant, while the bottom 
curve represents the minimal initial heating case, 
go = m c /m p . 

From Figure 6 it is evident that the electron and 
proton temperatures are indeed consistent under 



a planar model with some degree of initial elec- 
tron heating followed by Coulomb equilibration. 
The actual estimate of go will depend on which 
X-ray model is used to derive the electron tem- 
perature to compare to the proton temperature, 
as well as which region behind the shock is cho- 
sen. However, the relative electron temperatures 
and the ionization timescale differences across the 
three regions are simply fixed by the model as- 
sumptions. For a first estimate of go, we chose to 
use the go = m e /m p planar model to constrain T e 
at the outermost region (i.e., shortest ionization 
timescale, n c t). We compared this to the elec- 
tron temperature at that timescale implied by the 
proton temperature as a function of go- These es- 
timates are listed in Table 6. The errors include 
those due to T p , T c , and abundances summed in 
quadrature, which were all of approximately equal 
importance. As seen before from the comparison 
between T avCiS and T p , only X5, the aperture with 
the slowest shock speed, is consistent with go = 1. 
All others are consistent with go — m c /m p . The 
large errors in both T p and T c when added in 
quadrature lead to large errors in g ■ 

Given that our independent measurements have 
confirmed that the blast wave temperatures are 
consistent with a planar shock with some degree of 
initial electron-ion equilibration, it is now possible 
to constrain g by incorporating our knowledge of 
the proton temperature into the X-ray model fits. 
We expect this to greatly reduce the errors in our 
determination of g for two reasons. Firstly, T , 
T p , n c t and the abundances are all correlated vari- 
ables and hence summing their errors in quadra- 
ture over-predicts the actual error. Secondly, our 
three nested regions give us some constraint not 
only on the average electron temperature but also 
its spatial gradient. Allowing go to vary gives us 
access to different temperature evolutions, beyond 
the three extreme cases tested before. 

We adapted the planar model to consider a sin- 
gle T p and allow for different values of go from 
m c /m p to 1. For each aperture we built a grid of 
models for the best-fit T p as well as its la error 
limits. The only free parameters are then go and 
n c t . The average shock temperature and velocity 
are dependent on T p and go (see the above discus- 
sion of Figure 6). For any given g , the electron 
temperature varies with distance behind the shock 
as the downstream gas undergoes Coulomb colli- 
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sions. For consistency, the abundances and col- 
umn density values were fixed to the same values 
used in the previous fits. 

The fits to the variable go model are as good as 
the previous three. The resultant go values (listed 
in Table 7) are consistent with the previous cruder 
estimates but with tighter constraints. In many 
regions n c t and T p work together to give very sim- 
ilar g values within the la range on T p . One 
aperture, X4, is found to be of intermediate equi- 
libration, excluding both g — m c /m p and go = 1 
(at the lcr level). 

What do these g values imply for the rela- 
tionship between velocity and initial equilibration 
in collisionlcss shocks? In Figure 7 we compare 
our DEM L71 results on go as a function of V s 
with previous work on other remnants using opti- 
cal observations of the Ha line alone (Ghavamian 
et al. 2001). The results are consistent and in- 
dicate a decreasing level of equilibration with in- 
creasing shock speed. The method presented in 
the current paper is considerably different than the 
previous work modeling the Ha broad-to-narrow 
ratios. Ha emission is only produced extremely 
close to the shock front, so the go values modeled 
by this method truly reflect the initial equilibra- 
tion, modulo assumptions about the preshock ion- 
ization fraction. In contrast the X-ray emitting 
region is extended. We measure the electron tem- 
perature some distance downstream of the shock 
and then infer the initial equilibration, assuming 
that only Coulomb collisions have equilibrated the 
temperatures thereafter. That these two different 
methods lead to consistent go values as a function 
of shock velocity validates both methods. Fur- 
thermore it suggests that there are no significant 
additional heating mechanisms beyond Coulomb 
collisions that operate downstream of the immedi- 
ate shock zone. The real test, of course, will come 
when both methods can be applied to the same 
shock front. Unfortunately this was not possi- 
ble here, because the broad-to-narrow ratios were 
anomalously low, most likely due to precursor con- 
tamination of the narrow component of the Ha 
emission (see GRHW03 for more details). 

6. Summary & Concluding Remarks 

We have compared the post-shock electron 
and proton temperatures in DEM L71 to deter- 



mine the initial electron-ion equilibration, go = 
(T e /T p )o, over a range of shock speeds. First, we 
conducted independent measurements of T c from 
spatially resolved Chandra ACIS-S X-ray spectra 
of the outer rim, using three different noncquilib- 
rium ionization models for the post-shock temper- 
ature evolution. The first two models assumed a 
constant velocity shock, one with full initial equi- 
libration, go = 1, the other with only minimal 
initial equilibration, go = m c /m v , followed by 
post-shock Coulomb collisional equilibration. The 
third model was a fully equilibrated (g = 1) NEI 
Sedov solution. In all cases the ionization state 
was allowed to progress downstream. We com- 
pared the inferred average shock temperatures 
with the proton temperatures from the Ha line 
profiles reported in GRHW03. The slowest shock 
was consistent with full equilibration in both the 
constant velocity and Sedov models, but not with 
minimal heating. In contrast all the faster shocks 
were consistent with g — m c /m p but not g = l. 

Our independent measurements of T c and T p showed 
that all regions were consistent with a planar shock 
model, where the electron and proton tempera- 
tures started at some initial degree of equilibra- 
tion, g , and then slowly equilibrated downstream 
due to Coulomb collisions. Given this agreement, 
we constructed a planar X-ray model which in- 
corporated the measured value of T p and allowed 
for a variable go- With the additional constraints 
from the proton temperature we found one region, 
X4, with intermediate shock velocity for which 
both full and minimal equilibration could be ruled 
out, albeit only at la. Our results support the 
hypothesis from Ghavamian ct al. (2001), that 
there is an anti-correlation between shock speed 
and initial equilibration. 

The method presented here, of comparing the 
optically determined T p to the X-ray measured 
T c , is a powerful new technique for constrain- 
ing the initial temperature equilibration of SNR 
shock fronts. Its strength lies in the fact that 
neither the proton nor the electron temperature 
are strongly model dependent, as this work shows. 
On the other hand, DEM L71 was a particularly 
good remnant for this study. We do not expect 
cosmic-ray emission to contaminate our blast wave 
spectra because DEM L71 has the lowest radio 
flux of any of the 25 LMC SNRs in Mathewson 
ct al. (1983) and there is no evidence in the X- 



9 



rays for a high energy continuum from synchrotron 
emission (Hughes et al. 2003). Additionally, the 
clear separation between the ejecta and the blast 
wave allowed us to obtain blast wave spectra that 
are uncontaminated by SN ejecta. Other Balmer- 
dominated SNRs, such as Tycho or SN1006, may 
suffer more from an intermingling of ejecta into 
the blast wave zone, and complications to the fit 
from a cosmic ray synchrotron component. 

All the shocks we studied in DEM L71 are in 
the interesting velocity range where intermediate 
initial equilibrations appear to occur. However, 
due to the large errors on our current measure- 
ments of g , only one region was shown to be of 
intermediate equilibration, although two other re- 
gions were best fit by intermediate values of go- 
In the current fits these errors were dominated 
by the fitted error on the electron temperature, 
by the root-mean-squarc error on the abundances 
and column density, and to a lesser extent by the 
error in the proton temperature. We have been 
awarded a re-observation of DEM L71 with Chan- 
dra to quadruple the total exposure time. With 
this observation we will constrain the abundances 
of all elements with prominent lines separately for 
each aperture and obtain tighter limits on the elec- 
tron temperature. The increased statistics will 
also allow us to constrain the temperatures and 
timescales of each nested region separately. Hence 
we will be able to directly determine the temper- 
ature variation downstream of the shock, and ver- 
ify whether it follows the evolution given by our 
shock models. Additional optical spectroscopy of 
the non-radiative shock regions identified with the 
Fabry-Perot would be useful to reduce the uncer- 
tainties in the proton temperatures, and better 
sample the broad Ha component. 

We are grateful for the use of Leisa Townsley's 
CTI-correction algorithm and associated response 
matrices. We are pleased to acknowledge use- 
ful discussions with A. Sluis and T. B. Williams. 
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CER was supported by a NASA Graduate Student 
Researchers Program Fellowship. 
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Fig. 1. — Ha image of DEM L71. Rutgers Fabry- 
Perot (RFP) image of supernova remnant (SNR) 
DEM L71 at the wavelength of the narrow Ha 
line. Blast wave regions from which we extracted 
spectra are as marked. 

Fig. 2.— Chandra soft X-ray image of DEM L71 
(0.2—0.75 keV). X-ray spectra were extracted from 
three consecutive regions behind the shock at each 
aperture. 

Fig. 3. — Sedov temperature evolution. The Sedov 
solution for the ratio, T ave /T ave , s , of the average 
temperature of a parcel of gas to the current aver- 
age shock temperature function of the ratio, 
t/t c , of the evolving time to the current time. The 
evolution is plotted for the innermost parcel of ma- 
terial in each of the three regions behind the shock, 
labeled by their current fractional radii. Note that 
for the go — 1 case T c — T ave . 

Fig. 4. — Sedov ionization timescale evolution. 
The Sedov solution for the ratio of the evolving 
ionization timescale, n e t, to the product of the 
ambient electron density with the current age of 
the remnant. This is plotted as a function of t/t c 
for the innermost parcel of material in each of the 
three regions behind the shock, labeled by their 
current fractional radii. 

Fig. 5. — X-ray spectra of the blast wave. For each 
aperture we plot the Chandra X-ray spectra from 
three consecutive regions behind the shock, out- 
ermost to innermost, bottom to top. The model 
shown is a non-equilibrium ionization (NEI) pla- 
nar shock with minimal initial temperature equi- 
libration, g = m c /m p . 

Fig. 6. — Temperature evolution behind a planar 
shock. Comparison between the T p measured opti- 
cally and the T e measured in the X-rays for aper- 
ture XI. The data points represent T c for three 
consecutive regions behind the shock for the best- 
fit minimally equilibrated (go = m c /m p ), crosses, 
and fully equilibrated (go = 1), diamonds, pla- 
nar shocks. In both cases, the three data points 
were not measured independently but rather were 
fitted together for a single average temperature. 
The error bars include the la errors on T e as well 
as the RMS errors on the abundances and column 
density. The curves represent T c as a function of 
ionization timescale implied by the optically de- 



termined T p for various go'- 1.0, 0.8, 0.6, 0.4, 0.2, 
m e /m p , from solid to dashed. 

Fig. 7.— The degree of initial electron-proton 
equilibration. The ratio of the initial electron to 
proton temperatures, go, is plotted as a function of 
shock speed for the five regions in DEM L71 (solid 
error bars). These are consistent with the anti- 
correlation between go and shock speed found from 
modeling the broad-to-narrow flux ratios for the 
Balmer filaments in other remnants (dotted boxes, 
Ghavamian et al. 2001). The uncertainties in the 
equilibration for the DEM L71 regions include the 
la error bars on the best-fit go at the nominal 
T p , the range in go within the la error bars on 
T p , and the range in g within the RMS errors on 
the abundance and column density, all summed in 
quadrature. The range of velocities comes from 
the la range in proton temperatures at the best- 
fit go- For the other remnants the boxes represent 
the allowed range of velocities and equilibrations 
for which both the width of the broad line and the 
broad-to-narrow ratio could be modeled. 
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Table 1 



RFP apertures 


Chandra apertures 


3,4 


XI 


5, 6 


X2 


7, 8,9 


X3 


11, 12 


X4 


13, 14 


X5 



Table 2 

RFP Ha Results for Selected Regions in the DEM L 71 Blast Wave 





Ha profile fits 


(T e /T p ) 


= i 


(T e /T p ) = 


m c /m p 


Aperture 


Vfwhm (km s r ) 


Is /In 


V shock (km s _1 ) 


kT p (keV) 


V shock (km s _1 ) 


kT p (keV) 


1 

2 
3 
4 
5 


840±iiS 
98513° 

735l 8 °° 
450+4° 


u - oi -0.06 
U ' O4 -0.09 

n 4Q+ - 07 

66+ 08 
n 44+ 06 


1050±$ 
12401?,?° 
100511™ 
915l + |° 
555 7q 


1 2q+ u - 38 

i ' 8i -0.56 
i in+0.44 
i - i8 -0.32 

0.99l°J? 
0-36l°; ° 


olo -100 

96013* 

77C.+140 
' ' °-115 

710l 8 °° 
4301^° 


i oq+O- 38 

^y-o.so 

i 01 +0.90 
i - 8i -0.57 

1 in+0.46 

i - iS -0.32 

0.98l°-J? 
0.361°;^ 



Table 3 





Average A h and 


ABUNDANCES 




model 


10 2 ° atoms cm~ 2 




dex 


Ne 
dex 


Fe 
dex 


.go = 1, planar 
5o = P^nar 
go = I, Sedov 
combined 


5.8±1.8 
6.0±1.7 
7.1±2.3 
6.3±1.9 


8.391°;°^ 
8.40l°;i° 

8.411°;}° 

8.401°°^ 


7 R4+ - 11 
7 fifi+° 11 

< -DO-0.15 
7 7O+0.07 

7 7n +o.io 


6.86+°'°? 
6.861°;^ 

6.831°-°^ 
6.86+°'?! 
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Table 4 

Chandra X-ray Best-fit Model Parameters for Selected Regions in the DEM L71 Blast 

Wave 



Aperture 


NEI model 


kT ave s (keV) 


l0g(n c t fi na i) cm 3 s 


x 2 


9 

r-X 


i— 1 i— 1 i— 1 


.9o 
.90 = 
.90 


= 1, planar 

= ^ P lanar 
= 1, Sedov 


0.75lg;» 
0.95lg;g 


10.83t°;2o 
11.03lg;g 

10.34+g"™ 


129.45 
124.78 
143.93 


1.455 
1.402 
1.617 


2 
2 
2 


.90 
.9o = 
.90 


= 1, planar 
= planar 
= 1, Scdov 


U ' y ^-0.39 

1 2fi+°- 79 

1 - ZO -0.72 

o.74t° : « 


io-68±»:Sg 
10.97+^ 

10-23±g;Sl 


85.84 
79.76 
82.79 


1.431 
1.323 
1.380 


3 
3 
3 


.90 
.90 = 
.90 


= 1, planar 
= — , planar 
= 1, Sedov 


0.66±8;E 
70+ - 08 
n -54+0 05 


11.27±8;J« 

n 42 +0.06 


158.69 
144.36 
174.74 


1.430 
1.301 
1.574 


4 
4 
4 


9o 
.90 = 
.90 


= 1, planar 
= planar 
= 1, Sedov 


n 79 +U.2U 
u -' z -0.16 

0.94^ : I° 
0.66+^ 


10.37™ 


73.54 
78.15 
76.25 


0.994 
1.056 
1.030 


5 
5 
5 


.90 
.90 = 
.90 


= 1, planar 
= planar 
= 1, Sedov 


c ii «;;;;;;::;: 

u -^ z -0.03 
u - zo -0.02 


n.9« 
ii-93l8:S 
"■7218:5? 


164.85 
159.37 
152.21 


1.459 
1.410 
1.348 



Table 5 

Variation in kT c and n c t across the 3 Nested Regions Implied by each Shock Model 







outer region 


middle 


region 


inner region 


Aperture 


NEI model 




log(n e i) 


kT c 


log(n e t) 


kT c 


log(n e i) 






keV 


cm 3 s 


keV 


cm 3 s 


keV 


cm 3 s 


1 


go = 1, planar 


0.75 


10.35 


0.75 


10.65 


0.75 


10.83 


1 


90 = Planar 


0.60 


10.55 


0.73 


10.85 


0.80 


11.03 


1 


go = 1, Sedov 


0.72 


10.31 


0.96 


10.38 


1.25 


10.34 


2 


<7o = 1, planar 


0.94 


10.20 


0.94 


10.50 


0.94 


10.68 


2 


90 = Planar 


0.67 


10.49 


0.83 


10.79 


0.93 


10.97 


2 


<7o = 1, Sedov 


0.93 


10.20 


1.23 


10.27 


1.60 


10.23 


3 


go = 1, planar 


0.66 


10.79 


0.66 


11.09 


0.66 


11.27 


3 


•9o = T^S planar 


0.62 


10.94 


0.67 


11.24 


0.69 


11.42 


3 


<7o = 1, Sedov 


0.42 


11.11 


0.56 


11.18 


0.73 


11.14 


4 


50 = 1, planar 


0.72 


10.47 


0.72 


10.77 


0.72 


10.95 


4 


50 = Planar 


0.63 


10.63 


0.75 


10.93 


0.82 


11.11 


4 


<7o = 1, Sedov 


0.82 


10.34 


1.09 


10.41 


1.42 


10.37 


5 


go = 1, planar 


0.416 


11.42 


0.416 


11.72 


0.416 


11.90 


5 


50 = ^, Planar 


0.42 


11.45 


0.42 


11.75 


0.42 


11.93 


5 


go = 1, Sedov 


0.31 


11.69 


0.41 


11.76 


0.53 


11.71 
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Table 6 

Initial Electron-Ion Equilibration 



RFP Ha Chandra 

outermost region 



Aperture 




kT c 


log(n e t) 


.90 = (T c /T p )o 




(keV) 


(keV) 


cm 3 s 




1 

2 
3 
4 
5 


1 90+0.38 
1 81 +°- 95 

i - 8i -0.57 
i 1O+0.46 
i - iS -0.32 

0.98±g;» 
0.3618& 


o.eol^ 

0.67t°; 2 « 

o.62l8:S| 

u -°°-0.19 


10.55±»:ig 

io.94lg : gS 

10.631°^ 

ii.45l8:8S 


n 4l+ () - 31 

u - 4i -0.41 

o.33l°:i 

u - ZD -0.26 
u -°*-0.54 
1.0-0.13 



Table 7 
Variable g Model 



Aperture 


(km s- 1 ) 


fcTp (keV) 


50 


1 

2 
3 
4 
5 


occ + 120 
OOO-105 

980l 2 ?8 

77C + 140 
' ' °-110 

785l + | 


i 2Q+ - 38 

1 - zy -0.30 
i 01 +0.95 
i - 8i -0.57 
-i iD+0.46 
i - i8 -0.32 

o.98i8 : 2 2 ? 
o.36l8:S 


0.24+ - 12 

0.10+ 029 
0.01+003 

4fi+° 18 

u -*°-0.15 
1.0-0.16 
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